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^ ! ABSTRACT 

, ■ Aims. We carried out high-resolution spectroscopy and BV(I)c photometric monitoring of the two fastest late- type 

^ ■ rotators in the nearby /3 Pictoris moving group, HD 199143 (F7V) and CD-64°1208 (K7V). The motivation for this 

work is to investigate the rotation periods and photospheric spot patterns of these very young stars, with a longer term 

^-^, view to probing the evolution of rotation and magnetic activity during the early phases of main-sequence evolution. 

^^ i We also aim to derive information on key physical parameters, such as rotational velocity and rotation period. 

L^ ■ Methods. We applied maximum entropy (ME) and Tikhonov regularizing (TR) criteria to derive the surface spot map 

^i \ distributions of the optical modulation observed in HD 199143 (F7 V) and CD-64°1208 (K7 V). We also used cross- 

r^ , correlation techniques to determine stellar parameters such as radial velocities and rotational velocities. Lomb-Scargle 

O |. periodograms were used to obtain the rotational periods from differential magnitude time series. 

I ■ Results. We find periods and inclinations of 0.356 days and 21.5 deg for HD 199143, and 0.355 days and 50.1 deg for 

O I CD-64°1208. The spot maps of HD 199143 obtained from the ME and TR methods are very similar, although the 

^ , latter gives a smoother distribution of the filling factor. Maps obtained at two different epochs three weeks apart show 

^ . a remarkable increase in spot coverage amounting to ^ 7% of the surface of the photosphere over a time period of 

, ^, ' only ^ 20 days. The spot maps of CD— 64° 1208 from the two methods show good longitudinal agreement, whereas 

the latitude range of the spots is extended to cover the whole visible hemisphere in the TR map. The distributions 

^vj , obtained from the first light curve of HD 199143 show the presence of an extended and asymmetric active longitude 

^ . with the maximum filling factor at longitude ^ 325°. A secondary active longitude is present at ^ 100°. The spotted 

QQ • area distributions on CD— 64° 1208 show two active longitudes separated by about 180°, which is not unusual on such 

ryr\ ' very active stars. 

^ : Key words. Stars: rotation, Stars: late-type, Stars: individual: HD 199143, 00-64^208, Stars: activity. Stars: starspots 

O ■ 1. Introduction IZuckerman et all (|200lh and ISong et all (|2003|) report a 

^~~^ ^ ^ /iv T^ N 1 1. 1 f total of 35 known members of the BPMG, each with 

Moving groups (MGs are kinematically coherent groups of one or more characteristics indicative of extreme youth. 



X 



stars that probably share a common origin-the evapora- |Kaisler et al.l (12004) found a high lithium abundance in one 
tion of an open cluster, the remnants of a star formation ^f ^^^ gp^g members. The best-known members of the 



_ region, or a juxtaposition of several small star formation gp^g ^^^ early-type stars, whereas few studies have con- 

^ . bursts at different epochs in adjacent cells of the velocity ^entrated on the late-type stars representing three fourths 

>^ , field. Recent years have seen the discovery of several MGs ^f ^^^ ^^^^^ BP^g membership. The rotational proper- 

within lOOpc of the Sun (e.g.; TW Hydrae, Horologium- ^-^^ ^f th e late-type me mbers have been recently investi- 

Tucana and p Pictoris). Stars in these groups range be- ^^^^^ ^^ iMessina et al.l (l2010h as part of the RACE-OC 

tween millions and tens of million years m age, and their 13^0 iect 
study should tell us about the formation and evolution of 
sparse stellar associations. The stars themselves provide a 

rare glimpse of the very earliest main-sequence evolution We have carried out high resolution spectroscopy and 

and an opportunity to study magnetic activity at its ex- BV(I)c photometric monitoring of a selected sample of 

tremes. Young late- type fast rotating members of MGs offer young late- type stars (spectral type later than F2). The 

a unique laboratory for studying stellar dynamos, magnetic sample chosen for this work has been selected from previ- 

structures and coronal heating. ously established members of the BPMG, based on pho- 

The P Pictoris m oving group ( BPMG ) is one of the tometric and kinematic properties. The motivation for this 

youngest (~12 Myr, iKaisler et al.l |2QQ4|) and the clos - work is to investigate the rotation periods and photospheric 

est MGs to the Sun (~36pc, IZuckerman et al.l l2QQlf ). spot patterns of these very young stars, with a longer-term 
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view to probing the evolution of rotation and magnetic ac- 
tivity during the early phases of main-sequence evolution. 

Here, we present results for two members of our sam- 
ple, HD 199143 and CD-64°1208. These are the fastest 
rotators among the known late- type s tars in the BPMG . 
HD 199143 has a vsinz of 128 kms'^ (JTorres et all I2QQ6D 
which is the largest o f any known single solar-like star 
within 50 pc. van den Ancker et al] (|2QQQf ) proposed that 
HD 199143 has been spun up by accretion of material from 
a close T Tauri-like companion responsible for the emis- 
sion lines, the ultraviolet variability and the excess infrared 
emission. In the case of CD— 64° 1208, vsini > 100 km s~^ 
([Zuckerman et al.L l2001i) and its unusually strong Li for a 
10 Myr old star has be en postulated as possibl y being due 
to its very fast rotation (|Soderblom et al.lll993f ). Both stars 
represent a special opportunity to study the stellar atmo- 
spheric activity in young solar-like ultra- fast rotators. 

Further information on the targets and observations are 
given in Sect.[2]and Sect. [3] respectively. Stellar parameters, 
time series analysis and phasing of the light curves are de- 
scribed in Sect.m Sect. [5] and Sect. [6] respectively. The spot 
map distributions are shown in Sect. 4 refsec: spot map, while 
in Sect. [8] we present the results. The conclusions are given 
in Sect. H 



2. HD 199143 and CD-64°1208 

As members of the young P Pictoris moving group, the 
F7 dwarf HD 199143 {d = 47.7 pc) and the K8 dwarf 
CD-64°1208 {d = 29.2 pc) probe the immediate post- 
T Tauri phase, during which stars contract to the zero-age 
main sequence. HD 199143, is a wide (1^) binary system. It 
has the largest v sin i of any solar-like star known within 50 
pc {vs'mi = 128 km s-\ Torres et al., 2006). HD 199143 
forms a physical pair with the Li-rich late-type dwarf 
BD— 176128. It exhibits str ong chromospheric act i vity i n 
Ha and in the ultraviolet (|van den Ancker et al.L l2000f ). 
The metallicity of HD 199143, as estimate d from photo- 
metric indices in the survey of [Nordstrom etaTI ([2004), 
is [Fe/H]= —0.07. Its logarithmic X-ray luminosity (in 
erg s~^) as measured with RQSAT is l og L^ = 30.56, but 



Table 1. The log of the BV(I)c photometric observations of 
the targets from SAAO and SSO. 



log(Lx/i^boi) is -3.3 (IZuckerman et al". l200l| ) ^ whi ch puts 



the star below saturation (Jeffries & Maxted", '2005). If the 
star has a short period (<0. 5 days), it would be place d in 
the super saturation range ( Garcia- Alvarez et al.Ll2QQ8[ ). 

CD— 64° 1208 is a triple system consisting of a visual 
binary (whose primary has a K8 V spectral type as esti- 
mated from B V(RI)c photometry (Zuckerma n et all l200lf ) 
with a separat i on of 0.18 arcsec and /S.K = 2.3 mag 
([Chauvin et al.L |2010[ ) and a wide A7 V companion at a 
distance of 70 arcsec (HIP 92024). A spectrum obtained 
with the Double Beam Spec trograph (DBS) at the 2.3 m 
telescope in SSO ([Zuckerman et al.L 2001) is dominated by 
the K8 V component and shows strong Ha in emission 
(EW - 2.2 A), a strong Li 6708A feature (EW - 580 mA) 
and fast rotation. The unusually strong Li signature is pos- 
sibly due to its very fast rotation (e.g. see Soderblom et al., 
ll993D . 

In principle, the dilution of the optical flux of 
CD— 64° 1208 by its close visual companion may affect its 
light variations and have an impact on the spot modelling. 
If the companion is a late- type active star, say later than 
F5 V, and rotates with a period commensurable with that 
of the target star, its rotational modulation may produce a 



Object 


HJD 


Site 


No. of 




2453000.0+ 




datapoints 


CD-64°1208 


473.14-486.27 


SSO 


16 


CD-64°1208 


554.88-558.26 


SSO 


27 


CD-64°1208 


579.28-584.38 


SAAO 


57 


HD 199143 


555.02-558.32 


SSO 


27 


HD 199143 


578.34-584.52 


SAAO 


78 



spurious additional signal in phase with the light variation 
of CD— 64° 1208 that will introduce a systematic error on 
the derived distribution of its spotted area vs. longitude. 
An upper limit for this error can be estimated by assuming 
that the V — K colour of the visual companion is at least 
1.10 mag, c orresponding to a star of spectral type F5 V (cf. 
ICox L |2QQQI Sect. 7.5), and that its optical light modulation 
has an amplitude of 0.2 mag, i.e. comparable with that of 
CD— 64° 1208. Thus, the upper limit to the relative optical 
flux variation of the companion is AFy/Fy = 8.8 x 10~^, 
which implies a maximum systematic error in the evalua- 
tion of the total spotted area of the target of ~ 0.9 per cent 
of the star disc, on the hypothesis of completely dark spots. 
Since the companion is likely to be significantly cooler than 
F5 V, the actual error may be considerably smaller and may 
be safely neglected, given that the amplitude of the light 
variation of the target is ~ 0.2 mag in the Fpassband, cor- 
responding to a maximum spotted area of ~ 17 per cent of 
the star's disc (cf. Sect.[3|). 

3. Observations 

3.1. Photometry 

Photometric observations of the stars HD 199143 and 
CD-64°1208 were obtained with Cousins BV(I)c fil- 
ters during 2005 using telescopes at the South African 
Astronomical Observatory (SAAO) and the Siding Spring 
Observatory (SSO), Australia, and are summarised in 
Table [H SAAO observations were made using the Modular 
Photometer on the 0.5 m telescope at Sutherland. The 
observations were made with reference t o the Cousins 
UBV(RI)c standards in the E-regions (|Menzies et al.l 
[1989,). SSO observations employed the 0.6-m reflector and 
PMT detector. The measurements of targets, HD 199143 
and CD— 64° 1208, were alternated with two comparison 
stars. The light curves of HD 199143 and CD-64°1208 are 
shown in Fig. [T] and Fig. [2] respectively. 



3.2. Spectroscopy 

High resolution spectra of HD 199143 and CD-64°1208 
were obtained during August 2004 at SAAO, using the 
GIRAFFE flbre-fed echelle spectrograph attached to the 
1.9-m Radcliffe telescope. This instrument is a copy 
of the MUSICOS spectrograph currently in use on the 
Ber nard Lyot telescope at Pic du Midi Observatory, France 
(Ba udrand fc BohmlJl992f ). The resolving power is A/AA ~ 
42000, corresponding to a width of 2 pixels of the 1024 x 
1024 TEK CCD. A Th-Ar arc lamp was used for wave- 
length calibration. Arc spectra were taken at regular inter- 
vals to calibrate possible drifts. Flat-flelding was performed 
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Table 2. The log of spectroscopic observations of the targets and standard stars from SAAO. 



Object 


UT Start 


Exp. time 


No. of 
frames 


Comments 


CD-64°1208 


2004 Aug 03 


21:37 


1800 


1 


Target Star 


CD-64°1208 
HD 199143 
HD 4247 
HD 886 


2004 Aug 04 
2004 Aug 04 
2004 Aug 04 
2004 Aug 04 


20:51 
23:41 
02:12 
02:40 


1800 

1800 

600 

250 


1 
1 
4 
1 


Target Star 
Target Star 
FO V template 
B2 IV telluric std 


HD 142764 
HD 171391 


2004 Aug 06 
2004 Aug 06 


19:40 
20:32 


300 
1200 


2 
5 


K5 V template 

G8 HI radial velocity std 


HD 160032 
HD 199143 
HD 199143 


2004 Aug 07 
2004 Aug 07 
2004 Aug 07 


17:07 
19:36 
20:07 


600 
1200 
1200 


2 

1 
1 


F4V template 
Target Star 
Target Star 


HD 160032 


2004 Aug 08 


18:10 


500 


2 


F4V template 



Table 3. Stellar parameters 



Object 


B-V 


d 

(pc) 


vsini 
(km s-^) 


Vhel ± CTVhel 

(km s-^) 


U±au 
(km s-^) 


V±av 
(km s-^) 


W ±aw 
(km s-^) 


HD 199143 
CD-64°1208 


0.48 
1.2 


47.7 
29.2 


115.5±7.5 
121.3±15.3 


2.5it:2.9 

2.2±2.2 


-10.57±1.77 
-3.06zb2.06 


-15.83±1.03 
-11.10±1.33 


-9.45±0.90 
-15.63±1.80 



Phosed light curve of HD 199145 (P,„t = 0.5558 d) 






yi ^ m(. j2^a o ^ ^ 



Phosed light curve of CD-641208 (P,„t = 0.5545 d) 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 ] .: 

Phase 
Diomonds: dota with MJD < 53558.32; Triongles: 53578.34 < MJD < 53580.52; Asterisks: MJD > 53584.4 




Phase 

: doto with MJD < 53486.27; Triongles: 53554.87 < MJD < 53558.26; Asterisks: MJD > 53579.27 



Fig. 1. The V-band differential magnitude of HD 199143 Fig. 2. The V-band differential magnitude of CD-64°1208 
versus the phase as given by Eq. [H Different sym- versus the phase as given by Eq. [2j Different sym- 
bols refer to different epochs of observation as indicated bols refer to different observational epochs as indicated 
(MJD=HJD-24500000). (MJD=HJD-24500000). 



by taking the spectrum of a tungsten lamp and also by il- 
luminating the CCD with uniform light through a diffusing 
screen. The wavelength range was 4300-6680 A spread over 
52 orders. Table [2] gives the log of the observations. 

The spectra were reduced in a standard fashion, after 
which the continuum was normalized to unity. Spectral re- 
gions centred on some prominent lines of interest are shown 
in Fig. [31 together with the spectra of stellar "templates" of 
similar spectral type obtained during the same campaign. 
Highly broadened (FWHM > 250 km s~^) photospheric 
absorption lines (most prominently Ha, the Mg I triplet, 
the Na I doublet and H/3) are visible in the spectra of both 
HD 199143 and CD-64°1208. The same lines are present 



in the template stellar spectra, but are, as expected, much 
narrower. 

4. Stellar parameters 

4.1. Radial velocity and space motions 

Heliocentric radial velocities were determined using the 
cross-correlation technique. The spectra of each star were 
cross-correlated order by order, using the routine FXCOR 
in IRAF, against spectra of radial velocity standards 
of similar spectral types (see Table [2]). For each order, 
the velocity is derived from the position of the peak of 
the cross-correlation function (CCF) found by fitting 
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Fig. 3. High-resolution spectra of HD 199143 and CD— 64° 1208. For comparison, we also show the spectra of template 
stars with similar spectral types. 
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Fig. 4. The V-band differential magnitude of HD 199143 
versus time with overplotted a sinusoidal fit with the rota- 
tion period of 0.3558 d. Different symbols refer to different 
observatories, as labelled. 



a Gaussian to the top of the function. Radial velocity 
errors are computed by FXCOR based on the fitted peak 
height and the antisymmetric noise as described by 



iTonrv fc DavisI (|l979l ). The radial velocities calculated 
for each order are weighted by their errors, and a mean 
value is obtained for each observation. Orders including 
chromospheric features and prominent telluric lines were 
excluded when determining this mean velocity. Finally, a 
weighted mean radial velocity is determined using all the 
observations from all the observing runs. In Table [3] we 
Hst the average heliocentric radial velocity (T4iei) and its 
associated error (crvhei)- The final mean radial velocities 
we have obtained for HD 199143 and CD-64°1208 are 
Vr: = 2.5 and ^r = 2.2 km s~^, respectively. 



We have used these mean radial velocities, together with 
the spectroscopic parallaxes (d = 47.7 and 29 .2 pc) and 
prope r motions from the Tycho-2 Catalogue (|H0g et al.L 
120001 ). to calculate the Galactic space- velocity components 
(t/, y, W) in a right-handed coordinated system (positive 
in the directions of the Galactic center. Galactic rotation, 
an d the North Galactic Pole, respectively), as determined 
bv iMontes et al.l (|200lf ). The resultant values are given in 
Tabled The positions of HD 199143 and CD-64°1208 in 
the Galaxy velocity diagram confirm their membership of 
the BPMG. These value s are similar to those reported by 
(IZuckerman et al.L[2QQll ). 
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4.2. Rotational velocity 

To determine an accurate rotational velocity for HD 199143 
and CD -64° 1208, we again make use of the IRAF FXCOR 
cross-correlation routine. For each of the observing runs, 
the observed spectra of HD 199143 and CD-64°1208 were 
cross-correlated against the spectrum of the template star 
and the width (FWHM) of the cross-correlation function 
(CCF) determined. The calibration of this width to yield 
an estimate of v sin i is done by cross-correlating artificially 
broadened spectra of the template star with the original 
template star spectrum. Broadened template spectra are 
created for v^mi spanning the expected range of values 
by convolution with a theoretical rotational profile using 
the program STARMOD. The resultant relationship between 
V sin i and FWHM of the CCF was then fitted with a fourth- 
order polynomial. We have tested this method with stars 
of known rotational velocity and similar spectral type to 
HD 199143 and CD-64°1208. On average, we have ob- 
tained values with errors of about 5% compared with those 
from the literature. The uncertainties in the vsinz values 
obtained by this method have been calculated using the 
parameter R defined by Tonry & Davis (1979) as the ra- 
tio of the CCF height to the RMS antisymmetric compo- 
nent. This parameter is computed by the IRAF task FX- 
COR and provides a measure of t he signal-to-noise ratio of 
the CCF. iTonrv fc Da^dsl (|l979 t) show that errors in the 
FWHM of the CCF a re p ropor tional to (1 + i?)~\ and 
[Hartmann et al.l (|1986[ ) and lRhode et al.l (|2001) found that 
the quantity =bf sini(l + i?)~^ provides a good estimate for 
the 90% confidence level of a v sinz measurement. Thus, we 
have adopted ±'usini(l + i?)~^ as a reasonable estimate of 
the uncertainties in our 'L'sinz measurements. 

We have determined vsmi by this method for all the 
spectra available for HD 199143 and CD-64°1208. We 
found vsiiii values that match within the errors by us- 
ing spectra from different nights and at different epochs. 
We note that v sin i derived by this method is not sensitive 
to possible asymmetries of the CCF caused by starspots: 
such changes are mainly in the peak of the CCF, whereas 
the v^vai determination is based on the FWHM. The re- 
sulting error- weighted means for all the observing runs for 
HD 199143 and CD-64°1208 are 115.5±7.5 and 121.3±15.3 
km s~^, which are the values given in Table [3l Th e se va l- 
ues are very simila r to those given by iTorres et al.l (|2006f ) . 
IWeise et all (|2010l ) obtained a larger value of vsin i = 
161±12 km s-i for HD 199143. They cross-correlated the 
stellar spectra (3900-6800 A) with appropriate numeri- 
cal templates for th e respective stellar spectral type (see 
iBaranne et al.Lll979i) . 



5. Time series analysis: the search for rotation 
periods 

The spectral types and measured v sin i values confirm that 
HD 199143 and CD-64°1208 are fast rotators with periods 
Prot < 1 days. We therefore limited the search for peri- 
odicity in the V-band differential photometry to the range 
0.1 — 10.0 days to allow for possible modulations on shorter 
as well as longer periods . We applied the Lomb-Scargle peri- 
odogram (Scargle, 1982] for unevenly spaced data following 
the prescriptions of Home & Baliunas (1986). 

For HD 199143, the Lomb-Scargle periodogram shows 
the maximum power at a period of 0.3558 days with a false- 



alarm probability smaller than 10~^. The formal error on 
such a period is 0.004 days as derived from the FWHM of 
the corresponding periodogram peak, and is determined by 
the limited time extension of the data set (only 29.5 days). 
However, putting the V-band magnitude in phase with the 
0.3558-d period, we find that a period change as small as 
0.001 days gives rise to a significant increase in the dis- 
persion of the light curve points. Therefore, we prefer to 
assume the latter as the most probable error of the period 
determined. 

Visual inspection of the differential magnitude time se- 
ries (see Fig.H]) shows how good the agreement is between 
the data and the sinusoidal fit with the 0.3558-d rota- 
tion period. We have retrieved contemporary V-band ob- 
servations of HD 199 143 from t he ASAS archive (All Sky 
Automated Survey, (jPojmanslil . [l997) and these are plot- 
ted together with our own observations. Also, the ASAS 
data agree with the period we found. We have searched for 
the rotation period also in the ASAS time series. However, 
owing to the low precision of HD 199143, ASAS photome- 
try (~0.035 mag) and the low amplitude of the rotational 
modulation (~ 0.055 mag), our period search did not allow 
us to derive any significant periodicity. 

For CD— 64° 1208, the Lomb-Scargle periodogram ex- 
ibits two peaks of comparable power and significance level 
greater than 99%, one at P = 0.3545d and another at 
P = 0.5459d. However, only the shorter period, when com- 
bined with the estimated stellar radius {R = 1.1 i?©), 
gives an equatorial velocity consistent with the measured 
vsini = 121.3 km s~^. Moreover, our rotation period esti- 
mate is confirme d by the P = 0.3 45 ± 0.04 d value as de- 
termined by (Me ssina et al.l |2010[ ) making use of the ASAS 
photometry. In fact, in the case of CD— 64° 1208, although 
the ASAS photometry has a precision similar to that of 
HD 199143, the rotational modulation is a factor 3 larger 
(~ 0.16 mag), thus allowing them safely to establish the ro- 
tation period. A significant periodicity was found in 7 out 
of the 12 time segments in which the entire ASAS time se- 
ries was sectioned. The error in our derived period can be 
assumed to be of the order of 0.001 days because a larger 
uncertainty significantly increases the dispersion of the light 
curve points. On the other hand, the periodicity at 0.5459 d 
may be the result of aliasing and noise fluctuations affecting 
our time series. 



6. Phasing of the light curves 

The observations of HD 199143 are phased with the 
ephemeris: 



HJDn 



2453555.02423 + 0.3558 x E, 



(1) 



giving the epoch of the minimum light. The V-band light 
curve is shown in Fig. [TJ where different symbols refer to 
observations collected at different epochs in order to show 
how the shape of the light modulation changes vs. time. 
Specifically, three intervals can be identified during which 
the light curve remains stable, i.e. from HJD 2453555.02 
to 2453558.32 (diamonds), from 2453578.34 to 2453580.52 
(triangles) and from 2453584.43 to 2453584.52 (asterisks). 
The photometric points belonging to the first time interval 
sample the entire light curve, whereas those of the second 
interval show a gap between phases 0.2 and 0.5. A variation 
of the peak amplitude of the light curve is clearly evident 
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between the two time intervals, making it impossible to 
put the data together to improve the phase sampling of the 
light curve. The third data set does not show any significant 
variation and covers only the phases between 0.75 and 0.80: 
these points will not be considered for the subsequent spot 
modelling. The brightest differential F magnitude along the 
whole time series is -0.295 and was assumed as the unspot- 
ted magnitude for the spot modelling of HD 199143. An 
inspection of the ASAS magnitudes - which better sample 
the range of variability because of their 9 year baseline - 
confirms that our choice of unspotted value is the correct 
one. 

The observations of CD— 64° 1208 are phased with the 
ephemeris: 



HJD, 



2453473.20700 + 0.3545 x E, 



The V-band light curve is plotted in Fig. [2l with different 
symbols indicating data collected at different times. Also in 
this case, three time intervals can be identified during which 
the light modulation remained more or less stable, i.e. from 
2453473.14 to 2453486.27 (diamonds), from 2453554.87 to 
2453558.26 (triangles) and from 2453579.27 to 2453584.37 
(asterisks). The points of the first two data sets trace very 
similar light modulations, except for the phase interval be- 
tween 0.7 and 1.0. However, they were considered together 
in order to improve the phase sampling of the combined 
light curve that was fitted with our spot model. The third 
data set consists of points between phases 0.15 and 0.6 that 
are systematically brighter than the corresponding points 
of the first and the second data sets. Since the systematic 
difference is about 0.025 mag, the most likely explanation is 
a decrease in the area of the component of the spot pattern 
uniformly distributed over the star that made it brighter. 
A change of the area of the order of 2 — 3% of the star disc 
over a time interval of about 20 days is enough to account 
for such a variation and cannot be considered unusual for 
such a very active star. 

The brightest differential V magnitude throughout 
the whole time series is 0.829 and this was assumed 
as the unspotted magnitude for the spot modelling of 
CD— 64° 1208. As done in the previous case, we found that 
the correctness of this value is confirmed by the ASAS mag- 
nitudes timeseries. 



7. Spot map distributions 

7.1. Spot modelling technique 

The classic Doppler Imaging technique can provide us with 
a map of the stel l ar ph otosphere of an active star (e.g., 
IStrassmeierl, I2002L |2009[ ) but its applicability is strongly 
limited in the present case owing to the extreme rotational 
broadening of the photospheric spectral lines. Therefore, 
spot modelling based on wide-band light curve fitting ap- 
pears to be more appropriate for a preliminary characteri- 
zation of the photospheric activity of our stars. 

The reconstruction of the brightness distribution over 
the surface of a star by modelling the rotational modula- 
tion of wide-band optical fluxes is an ill-posed problem. The 
modulation provides information only on the variation of 
the projected spot area versus the rotation phase, i.e., the 
stellar longitude. In principle, it is possible to determine a 
unique spot distribution by minimizing the y^ of the resid- 
uals between the observed and the model light curve, but 



such an approach is unsatisfactory because the solution is 
unstable, i.e. small variations in the input data lead to large 
changes in the spot map. This is due to the role played by 
the noise in the y^ minimization, implying that most of 
the structure appearing on the spot map actually comes 
from the overfitting of the light fluctuations produced by 
intrinsic stellar variability or measurement errors. 

It is possible to overcome the uniqueness and stability 
problems by introducing a regularizing function into the 
solution process. This corresponds to the a priori assump- 
tion of some specific statistical properties for the spot map, 
which allows us to select one stable map among the poten- 
tially infinite maps that can fit a given light curve. The two 
most used a prio ri assumptions a r e the maximum entropy 



(hereinafter ME; 'Gull & Skillini, Il984l: IVogt et al. 



V / and the Tikhonov (hereinafter TR; iPiskunov et al. 



1987 ) 



19901) 



regularizations. These are particularly well-suited when the 
spot map consists of an array that specifies the spot cov- 
ering factor in each surface element of the model star. The 
covering factor / gives the specific intensity / of each sur- 
face element according to the definition: / = (1 — /)/u + /^s5 
where I^ is the specific intensity of the unspotted photo- 
sphere of the surface element, /g that of the spotted photo- 
sphere, and < / < 1. 

A detailed description of the application of the ME and 
TR regul arizations to spot modelling problems was pre- 
sented by iLanza et al.l (|1998[ ) to whom we refer the reader 
for a detailed discussion of our approach. Here we only re- 
call that the ME spot maps are computed by a constrained 
minimization of a functional Qme, 



mE 



X - Ame^*, 



(3) 



whereas the TR maps are computed by the constrained 
minimization of 



?TR = X + XtrTR. 



(4) 



Here, x^ is the normalized sum of the squared residuals 
between the observed and the synthesized light curves, S 
is the entropy of the spot map, TR is the Tikhonov func- 
tional of the spot map, and Ame > and Atr > are the 
Lagrange multipliers for the ME and the TR regulariza- 
tions, respectively. The explicit expressions for S and TR, 
and the procedu re of evaluation of th e Lagrange multipliers 
are described by Lanza e t al.l (|l998f ) . 

The synthesized light curve for a given distribution of 
the covering factor is computed by adopting a spheroidal 
geometry for the star, i.e. it is represented as an ellipsoid 
of semi-axes a = b and c < a, whose flattening is due to the 
centrifugal potential produced by the rapid rotation. We 
compute the flattening in the Roche approximation as 



1 



GM 



(5) 



where c is the polar radius, a the equatorial radius, Q the 
angular velocity of rotation, M the mass of the star and G 
the gravitation constant. The effective temperature T on 
the photosph ere of a distorted star d epends on the local 
gravity (e.g., iKallrath fc Milona , Il999 ). For a star having 
an outer convective envelope, a good approximation is given 

by 



r = r„ 



1 + 0.08 X 



9p 



(6) 



D. Garcia- Alvarez et al.: Starspots on the fastest rotators in the /3 Pic moving group 



where Tp is the effective temperature at the pole, g is the 
local surface gravity and g^ is the gravity at the pole. By 
considering such a gravity darkening effect, we compute the 
distribution of the effective temperature on the surface of 
the star and adopt Kurucz, (;2QQQ,) models to compute the 
brightness emerging from each surface element. Adding the 
effect of the limb-darkening, treated in the linear approxi- 
mation, we derive the unperturbed brightness lu over the 
stellar photosphere. The brightness of the spotted photo- 
sphere is given by: /g = Cg/u, where the contrast C^ de- 
pends on the spot effective temperature and is assumed 
constant over the photosphere of the star. In order to com- 
pute the emerging flux with a relative precision of the order 
of 10~^, the stellar photosphere is subdivided into squared 
elements of side Sq = 1°, the fluxes of which are summed 
to compute the total flux from the stellar disc. 

Absolute properties of the spots cannot be extracted 
from single-band data because systematic and other errors 
arise from the unknown unspotted light level of the star and 
the assumption of single-temperature spots. Specifically, 
the brightest magnitude observed in such active stars is 
probably fainter than the true unspotted magnitude be- 
cause spots may always be present on their photospheres. 
However, our model can be applied to derive a map of 
the component of the spot pattern un evenly distributed 
vs. longitude (for example, see Lanza etall , l2QQ6l for 
a detailed discussion and a comparison with the capabil- 
ities of Doppler Imaging maps). The uniformly distributed 
component of the spot pattern cannot be derived because 
we assume that the unspotted magnitude is that corre- 
sponding to the brighter hemisphere of the star. Only when 
a long-term sequence of seasonal light curves is available 
can we extract information on the variation of the uni- 
formly distr i buted spot pattern from season to season (cf. 
iLanzaet alll2QQ2[ ). 

7.2. Model parameters 

In our spot maps, we consider squared surface map elements 
of side s = 18°. while fluxes are always computed with the 
finer subdivision Sq = 1°. The stellar parameters adopted 
for the modelHng of the two objects are listed in Table H] 
together with the reference from which they were taken. 
The spot temperatures are estimate d from the B — V and 
V — I colour variations accordin g to iMessina et al.l (|2QQ6al ) 
and lMessina fc GuinanI (|2QQ6b[ ). 

Specifically, we model the amplitude of the V magni- 
tude versus B — V and V — I colour variations induced 
by the cool/hot spots on the stellar photosphere for a 
grid of spot tem peratures and filling factors using the 
Dorren app roach (DorrenlJ 19871 ) and the Hauschildt et al. 
(IHauschildt et al. I Il999[ ) atmosphere models. A x^ mini- 
mization is used to select the solutions that best fit the 
observed amplitudes. As expected, the solutions are not 
unique and the higher the spot temperature, the higher 
their filling factor values. However, the solutions tend to 
cluster around a finite range of values for the spot temper- 
ature. Such a finite range of temperature values is adopted 
as the uncertainty of our spot temperature determination, 
which is about ±150 K for both stars under analysis. 
Assuming that the observed magnitude and colour varia- 
tions arise from both cool and hot spots, what we infer from 
our modelling is an average value for the temperature inho- 
mogeneities. In the case of HD 199143 the spot temperature 



Table 4. Stellar parameters of HD 199143 
CD— 64° 1208 adopted for our spot modelling. 



and 





HD 199143 


CD-64°1208 


Ref. 


SpTy 


F8V 


K8V 


1 


M (Mq) 


1.50 


0.80 


2 


a{RQ) 


2.20 


1.10 


2 


Prot (d) 


0.3558 


0.3545 


3 


i (deg) 


21.5 


50.1 


3 


c/a 


0.622 


0.911 


3 


Teff (K) 


6310 


4260 


1 


Ts{K) 


4700 


3700 


3 


\ogg (cm s ^) 


3.9 


4.3 


2 


Uv 


0.667 


0.800 


3 


a 


0.234 


0.315 


3 



References: 1. Allende Prieto & Lambert (1999); 2. de La Reza 
& Pinzon (2004); 3. This work. 



inferred from either 5 — V or V — / is in agreement within 
the computed uncertainty. For CD— 64° 1208, the spot tem- 
perature derived from the B — V colour is about 200 K 
warmer than that derived from the V — I colour, probably 
because of the excess flux of faculae mainly affecting the 
B-band (jMessina et al.Ll2QQ6ai ). 

Therefore, the V — I colour probably gives a better de- 
termination of the cool spot temperature; we estimate an 
uncertainty of ±150 K for both stars. It is important to note 
that cool spots dominate the flux variation of both stars in 
the V-band whereas the possible facular component can be 
regarded as a second order effect in such extremely active 
stars (cf. th e case of HR 1099 as analysed and discussed by 
IVogt et aD , 11999). Therefore, the temperature of the cool 
spots is the only one appropriate for our modelling of the 
optical flux rotational modulation. From the effective tem- 
peratures of the spotted and the unperturbed photospheres, 
the spot contrast Cg = Is/Iu at the disc centre is evaluated 
by means of Kurucz models, adopting the same mean grav- 
ity (see the values listed in Table HJ. 

The linear limb-darkening coefficients in the V-band, 
iiy , are assumed to be the same for the spotted and unspot- 
ted photospheres, which does not introduce any significant 
systematic error given that the difference between the two 
coefficients is less than 15% , and that the spots occupy a 
maximum fraction of the order of 10% of the stellar disc. 

The inclination i of the stellar rotation axis with respect 
to the line of sight is derived from the observed vsini, the 
estimated stellar equatorial radius a and the rotation pe- 
riod Prot found in Sect. [H Given the errors in the mea- 
surements of vsini, the error in the inclination is of the 
order of 10° — 15°. Nevertheless, changing the inclination 
does not significantly affect the distribution of the relative 
spotted area vs. longitude that is the primary result of our 
modelling. 

The rapid rotation of the two stars produces a remark- 
able centrifugal fiattening that makes it necessary to in- 
clude gravity darkening and ellipsoidicity in the stellar 
model. Any systematic error in the estimation of such 
effects modifies the unperturbed fiux level of the star. 
However, to first order there is no difference in the spotted 
area between the most and the least spotted hemispheres, 
and the component of the spotted area unevenly distributed 
in longitude can be safely mapped. 
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Fig. 5. Upper panel: The first light curve of HD 199143 
(filled dots) with the corresponding synthetic light curve 
from the ME best fit spot distribution (solid line). The flux 
is given in relative units, assuming as the reference value 
(1.0) the brightest observation of the star (see Sect. 2). 
Middle panel: the residuals between the observations and 
the model in relative units. Lower panel: The ME spot map 
of the star in orthographic projection, showing the disc of 
the star as it would be seen by an observer on the Earth at 
the labelled rotation phases (see the text). 

8. Results 

The results of our spot modelling are summarized in 
Table [5l The reduced x^ is deflned as: 



Fig. 6. As Fig. [5l but for the T-regularized best fit. 



X 



M 

M ^ 



{0,-F,f 



ij) 



where Oi is the observed flux and Fi the computed flux at 
the z-th rotational phase. The standard deviation used to 
compute the best flt for each light curve has been derived 
a posteriori by considering the standard deviation of the 
differences between the observations and the best flt flux 
values, both normalized to the maximum observed flux. It 
is comparable with the a priori standard deviation expected 
from the accuracy of the photometry for the flrst light curve 
of HD 199143 and the light curve of CD-64°1208 (i.e. 0.005 
mag), whereas it is about 2 times smaller for the second 
light curve of HD 199143. 

The best flt of the flrst light curve of HD 199143 is shown 
in Fig. O together with the ME spot map in orthographic 
projection as seen by an observer looking at the star from 
the Earth. The north pole is always in view because of 
the small inclination of the rotation axis of the star. The 
reference frame is chosen so that the longitude of the central 
meridian of the stellar disc at rotation phase (j) (with < 
(^ < 1) is 360° X (j). The T-regularized best flt and the 
corresponding spot maps are shown in Fig. [6l The ME- 
and the T-regularized best flts of the second light curve of 
HD 199143 and of the light curve of CD-64°1208 are shown 
in Figs [TlfTOl respectively, together with the corresponding 
spot maps. 

The spot maps of HD 199143 obtained with the ME 
and the T regular izat ions are remarkably similar, although 
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Fig. 11. The distribution of the spotted area vs. longitude 
for the ME (solid line) and the TR (dotted line) models 
of the flrst light curve of HD 199143. The area unit is the 
surface of the star. The distributions have been plotted be- 
tween —60° and 420° to better display the structure of the 
spot pattern; the longitude interval 0° — 360° is marked by 
two vertical dashed lines. 



the TR criterion gives a smoother distribution of the fllling 
factor, as expected on the basis of its a priori assumptions 
that select the smoothest spot distribution that flts the 
data. The latitude of the spots is conflned below ~ 60° in 
the flrst map and below ~ 30° in the second map. This is 
due to the low inclination angle of the rotation axis that 
makes it necessary to put spots away from the north pole 
(which is always in view) to produce a sizeable magnitude 
variation vs. stellar rotation phase. The second light curve 
has a larger amplitude than the flrst one, making it neces- 
sary to increase the fllling factor and put the spots at lower 
latitudes to reproduce the observed amplitude. The lack of 
spotted area in the ME map between longitudes 100° and 
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Table 5. Results of the spot modelling of the light curves of HD 199143 and CD-64°1208. 



Star 



initial HJD final HJD 



M 



Ame Atr 



Xme 



Xtr 



^ME 



Atr 



HD 199143 
HD 199143 
CD-64°1208 



2453555.02 
2453578.34 
2453473.14 



2453558.32 
2453580.32 
2453558.26 



27 
63 
43 



0.00620 
0.00246 
0.00683 



1.0 

0.15 

2.5 



40.0 

6.0 

50.0 



1.1222 
1.0366 
0.8542 



1.0967 
1.0506 
0.8529 



0.0200 0.0678 
0.0910 0.2392 
0.0720 0.1620 



From left to the right, the name of the star, the initial and final HJD of the light curve, the number of photometric points in the light curve M, their 
standard deviation a, the Lagrangian multiplier of the ME and TR optimizations, the corresponding reduced % and the total spotted area in units of 
the entire photosphere of the star, respectively. 
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Fig. 12. As Fig. [TTJ but for the second light curve of 
CD-64°1208. 
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Fig. 13. As Fig.ini but for the light curve of CD-64°1208. 



200° is due to the lack of data in the second light curve 
to constrain the spot pattern in the corresponding phase 
range. We note the remarkable increase in the total area of 
the spot component unevenly distributed in longitude re- 



quired to fit the light curve amplitude variation. The ME 
solutions give an increment of ~ 7% of the entire photo- 
sphere on a time scale of only ~ 20 days. 

The spot maps of CD-64°1208 obtained with the ME 
and the T regularizing criteria show good longitudinal 
agreement whereas the latitude range of the spots is ex- 
tended to cover the whole visible hemisphere in the TR 
map. This is due to the lack of latitudinal constraints on 
the spot location due to the higher inclination of the stellar 
rotation axis. Therefore, in the case of CD— 64° 1208 what 
we can retrieve from the light curve modelling is basically 
the distribution of the non-uniform component of the spot- 
ted area vs. longitude. 

The distributions of the spotted area vs. longitude are 
shown in Figs [TT1IT21 for the light curves of HD 199143 and 
in Fig. [13] for the light curve of CD— 64°1208, respectively. 
The distributions obtained with the ME and TR methods 
are always very similar although there is a systematic dif- 
ference between the ME and the TR areas, the latter being 
systematically greater than the former. This is due to the 
different assumptions of the two regularization methods, 
with ME trying to decrease the spotted area as much as 
possible while TR tries to make the spot distribution as 
uniform as possible. Therefore, the TR maps have spots 
in the invisible part of the stellar surface that systemati- 
cally increase the total spotted area in each longitude bin. 
Of course, the spotted area projected onto the visible disc, 
which is responsible for the flux modulation, is the same 
for both solutions. 

The distributions obtained from the first light curve of 
HD 199143 show the presence of an extended and asym- 
metric active longitude with the maximum filling factor at 
longitude ~ 325°. A secondary active longitude is present 
at ~ 100°. It appears to have grown significantly in the 
second light curve maps while the primary active longi- 
tude has migrated slightly towards decreasing longitudes. 
Unfortunately, the incompleteness of the second light curve 
makes such results on the spot pattern evolution quite un- 
certain. 

The spot distribution derived for HD 199143 and 
CD— 64° 1208 can be compared with those obtained for the 
two fast-rotating main-sequence stars AB Dor and LO Peg. 
AB Dor is a KO dwarf with a rotation period of 0.5148 days 
that h as been extensi v ely st udied through spot modelling 
(e.g. 'Jarvinen e t al.L l2005f ) and Doppler Imaging (e.g. 



Jeffers et al., 2007|, and references therein). Spot modelling 
based on long-term seasonal photometry generally show one 
or two active longitudes that change from season to season. 
In Doppler images, it shows an asymmetric spotted polar 
cap of varying extension in addition to several low and in- 
termediate latitude spots. 

LO Peg is a K5-K7 dwarf with a rotation period of 
only 0.4236 days that has been Doppler imaged by, e.g.. 
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iBarnes et aD (|2QQ5[ ) and iPiluso et all (|2QQ8[ ) showing an 
asymmetric polar cap with several appendages reaching 
10° — 20° of latitude. Therefore, the spot distributions ob- 
tained from our limited photometric datasets are in general 
agreement with those obtained for other rapidly rotating 
young late- type stars. The presence of a more or less asym- 
metric polar cap can be established only through Doppler 
imaging because it does not contribute significantly to the 
optical flux modulation that is dominated by intermedi- 
ate and low-latitude spots. Of course, Doppler images have 
a greater spatial resolution than the maps obtained from 
the inversion of optical photometry and show that spots as 
small as a few degrees of longitudes are present on AB Dor 
and LO Peg. Therefore, our active longitudes are likely to 
be clusters of much smaller individual spots. 



asymmetric active longitude with the maximum filling 
factor at longitude ~ 325°. A secondary active longi- 
tude is present at ~ 100°. It appears to have grown 
remarkably in the second light curve maps, while the 
primary active longitude has migrated slightly towards 
decreasing longitudes. The spotted area distributions on 
CD— 64° 1208 show two active longitudes separated by 
about 180° which is not unusual on such very active 
stars. 
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9. Conclusions 

We have analysed photometric and high resolution spec- 
troscopic observations, made during the period 2004-2005, 
to investigate the surface activity of the two fastest rotat- 
ing late-type stars known in the P Pictoris young moving 
group. In addition to surface spot map distributions, these 
observations have yielded information on the key physical 
parameters of rotational velocity, rotation period, inclina- 
tions angle and radial velocity. 

1. Radial velocity standard stars have been used to calcu- 
late accurate heliocentric radial velocities for each ob- 
servation by using the cross-correlation technique. The 
mean velocities for HD 199143 and CD-64°1208 are 2.2 
± 2.9 km s~^ and 2.5 ± 2.2 km 5~^, respectively. These 
values were used to calculate the Galactic space- velocity 
components {U^V^W) and confirmed their membership 
of the BPMG. 

2. From cross-correlation analysis of our spectra, projected 
rotational velocities, vsinz, of 115.5±7.5 km s~^ and 
121.3±15.3 km s~^ have been measured for HD 199143 
and CD-64°1208, respectively. 

3. Rotation periods of 0.356±0.004 days for HD 199143 
and 0.355±0.040 days for CD-64°1208 have been de- 
rived from the time series of V-band differential magni- 
tudes. 

4. The inclination, i, of the stellar rotation axis with re- 
spect to the line of sight has been derived from the ob- 
served vsini, the estimated stellar equatorial radius a, 
and the rotation period Prot • We find inchnation angles 
of 21?5 and 50° with an uncertainty of 10° — 15° for 
HD 199143 and CD-64°1208, respectively. 

5. The spot maps of HD 199143 obtained using the TR and 
ME criteria are remarkably similar; the TR criterion 
gives a smoother distribution of the filling factor, as 
expected. The latitude of the spots is confined below 
^ 60° in the first map and below ~ 30° in the second 
map. Light curve amplitude variations over the different 
epochs of observation require an increase of ~ 7% in the 
spot filling factor, unevenly distributed in longitude, on 
a time scale of only ^ 20 days. 

6. The spot maps of CD— 64° 1208 show good longitudinal 
agreement whereas the latitude range of the spots is 
extended to cover the whole visible hemisphere in the 
TR map. 

7. The distributions obtained from the first light curve 
of HD 199143 show the presence of an extended and 
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Fig. 7. As Fig. O but for the second Hght curve of HD 199143. 
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Fig. 8. As Fig. ^ but for the second hght curve of HD 199143. 
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CD-641208 - (HJD 2453473.14 - 2453558.26) - MEM map 
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Fig. 9. As Fig. [51 but for the Hght curve of CD-64°1208. 
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Fig. 10. As Fig. El but for the hght curve of CD-64°1208. 



